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Abstract 

o 

■ The theory of QED corrections to hyperfine structure in light hydrogenic 

. atoms and ions has recently advanced to the point that the uncertainty of these 

"p^ | corrections is much smaller than 1 part per million (ppm), while the experi- 

ments are even more accurate. The difference of the experimental results and 
the corresponding QED theory is due to nuclear effects, which are primarily 
the result of the finite nuclear charge and magnetization distributions. This 
difference varies from tens to hundreds of ppm. We have calculated the domi- 
nant nuclear component of the Is hyperfine interval for deuterium, tritium and 
singly ionized helium, using a unified approach with modern second-generation 
potentials. The calculated nuclear corrections are within 3% of the experi- 
mental values for deuterium and tritium, but are roughly 20% discrepant for 
helium. The nuclear corrections for the trinucleon systems can be qualitatively 
understood by invoking SU(4) symmetry. 
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1 Introduction 



Until very recently hyperfine splittings in light hydrogenic atoms were the most pre- 
cisely measured atomic transitions. Theoretical predictions based on QED are less 
accurate, but have improved considerably in recent years. Non-recoil and non-nuclear 
contributions [1, 2] are known through order a 3 Ep, where Ep is the Fermi hyperfine 
energy (viz., the leading-order contribution) and a is the fine-structure constant. Be- 
cause the hadronic scales for recoil and certain types of nuclear corrections are the 
same, recoil corrections are treated on the same footing as nuclear corrections [1], and 
we will call both types "nuclear corrections." Uncalculated QED terms of order a 4 E F 
in light atoms are almost certainly smaller than .1 ppm., while the experimental er- 
rors are smaller still. This provides us with the unprecedented opportunity to study 
nuclear effects in the hyperfine structure (hfs) of light hydrogenic atoms, which range 
in size from tens to hundreds of ppm. We will restrict ourselves to hydrogenic s-states, 
because these states maximize nuclear effects. 

Table 1: Difference between hyperfine experiments and QED hyperfine calculations 
for the nth s-state of light hydrogenic atoms times n 3 , expressed as parts per million 
of the Fermi energy. This difference is interpreted as nuclear contributions to the hy- 
perfine splitting [2] . A negative entry indicates that the theoretical prediction without 
nuclear corrections is too large. 



n 3 (E^~E^)/E F (ppm) 



State 


H 2 H 3 H 3 He+ 


Is 


-33 138 -38 -212 


2s 


-33 137 - - -211 



Table 1 is an updated version of the corresponding table in Ref. [2] . Because nuclear 
effects have a much shorter range than atomic scales, we expect the splittings in the 
nth s-state to be proportional to |0 n (O)| 2 ~ 1/n 3 , where 4> n (r) is the non-relativistic 
wave function of the electron. Forming the fractional differences (in parts per million) 
between E^ and E^ D (times n 3 ) leads to the tabulated results. As we stated above, 
these large differences reflect neither experimental errors nor uncertainties in the QED 
calculations, but rather the large nuclear contributions to hfs. 

In order to perform a tractable calculation it is necessary to restrict the scope, 
while at the same time incorporating the dominant physics. To accomplish this we 
borrow a technique from chiral perturbation theory (%PT, the effective field theory 
for nuclei based on QCD). This technique, called power counting, is the organizing 
principle of xPT and allows one to perform systematic expansions [3] in powers of a 
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small parameter, (Q/A), where Q is a typical nuclear momentum scale that can be 
taken to be roughly the pion mass (m^ ~ 140 MeV), and A is the large-mass QCD 
scale (~ 1 GeV) typical of QCD bound states such as the nucleon, heavy mesons, 
nucleon resonances, etc. We also note that 1/Q specifies a typical correlation length 
(and a reasonable nearest-neighbor distance) in light nuclei (~ 1.4 fm)[4]. This ex- 
pansion in powers of (Q/A ~ 0.1 — 0.15) should converge moderately well. In this 
work we restrict ourselves to leading order in this expansion, and this restriction elim- 
inates nuclear corrections of relativistic order, which are subleading and exceptionally 
complicated because of the complexity of the nuclear force [5, 6]. 

In processes that involve virtual excitation of intermediate nuclear states (each 
with its own energy, E N , relative to the ground-state energy, E ) the excitation 
energy (oun — -EW — -^o) is of order Q 2 / A and typically is a correction to the leading 
order [4]. Consistency therefore demands that we drop such terms. The nuclear recoil 
energy scales like Q 2 /M, where M is the nucleon mass, and can also be dropped. 
These are very considerable simplifications in constructing the nuclear contribution 
to hfs, which we present in the next section. 



2 Nuclear Contributions to Hyperfine Structure 

The hyperfine interactions that interest us are simple (effective) couplings of the 
electron spin to the nuclear (ground-state) spin: a ■ S, where a is the electron (Pauli) 
spin operator and S is the nuclear spin (total angular momentum) operator. Other 
couplings of the electron spin are possible and either generate no hyperfine splitting, 
none in s-states, or higher-order (in a) contributions. 

The lowest-order (Fermi) hyperfine interaction is generated by the interaction of 
the electron current with the magnetic dipole part (determined by pi) of the nuclear 
current. A simple calculation gives the well-known result [1] 

_ 47ra/i 7V |0 n (O)| 2 3 ■ S 
hi-p = — 



3m 



e 



where (SS\fj, z \SS) = /ijv defines the nuclear magnetic moment and m e is the electron 
mass. The factor of (a-S)/S leads to a hyperfine splitting proportional to (2S+\)/S. 
All additional contributions will be measured as a fraction of this energy. 

Naively calculating the higher-order (in a) corrections using only the nuclear mag- 
netization distribution will fail. The atomic wave function is modified by the nuclear 
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(a) (b) (c) 



Figure 1: Nuclear Compton amplitude with direct (a), crossed (b), and seagull (c) 
contributions illustrated. Single lines represent an electron, double lines a nucleus, 
and shaded double lines depict a nuclear Green's function containing a sum over 
nuclear states. The seagull vertex maintains gauge invariance. The four-momentum 
running through the loop is g M . 

charge distribution in precisely the same region that the magnetization distribution is 
nonvanishing[7, 8]. It is therefore necessary to incorporate the complete set of second- 
order (in perturbation theory in a) processes shown in Fig (1), which comprise the 
nuclear Compton amplitude coupled to the electron Compton amplitude. Only the 
forward-scattering part of this amplitude is required for the 0(a) corrections to E-p. 
The resulting energy shift is then given by 

AE = .(4™) |0 n (O)| J — ^ - _ — - - ) , (2) 

where t^ u is the electron Compton amplitude and T^ v is the nuclear Compton ampli- 
tude, which is required to be gauge invariant. The lepton tensor can be decomposed 
into an irreducible spinor basis, and we can ignore odd matrices and spin-independent 
components. 

We also ignore (for now) terms that couple two currents together. It is easy 
to show that since the nuclear current scales as 1/A (the conventional components 
of the current have explicit factors of 1/M), two of them should scale as 1/A 2 and 
generate higher-order (in 1/A) terms. This leaves a single dominant term representing 
a charge-current correlation 

a p U ,2 U m2 [ d 4 q (ffx qr[T m0 (q,-q)-T 0m (q,-q)] 

AE = (Ana) |0 n (O)| J — ^ - _ — - - } . (3) 

The nuclear seagull terms B 0m (q, —q) and B m0 (q, —q) that are contained as part of 
T^ u in Fig. (lc) are of relativistic order[5] (~ 1/A 2 ) and can be dropped. Although 
the seagull terms B mn (q, — q) are of non-relativistic order, they do not contribute to 
hfs because of crossing symmetry. The explicit form for the remaining term in T m0 
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(suppressing the nuclear ground-state expectation value), which involves a complete 
sum over intermediate states, N, is given by 

r „„ (9i , = E j^tMmm + crossed tem , (4) 

where J(q) and p(q) are the Fourier transforms of the nuclear current and charge 
operators. The crossed term has the operator order reversed and q — > — q . 

The limits uj n — > and m e — > (both scales are small compared to the nuclear-size 
scale, l/Q) greatly simplify the calculation, and this leads to 

AE = i(4na) |0 n (O)| J > ( 5 ) 

which is infrared divergent. Using 

J m (- q)= J d 3 y J m (y) exp {-i q ■ y) , (5) 

and 

p(q) = / d 3 x p(x) exp (i g • x) , (6) 
together with z = x — y, and a lower-limit (infrared) g— cutoff, e, we find 



A£ = -8a 2 |0 n (O)| 2 J d*x J d*y {p(x), 3 ■ J(y)} x V,(^ - z - + ^-) , (7) 

where there is an implied (nuclear and atomic) expectation value. The constant term 
does not contribute because of the derivative. The second term is only logarithmically 
divergent when Siegert's theorem is applied[9], can be shown to vanish in several limits 
(including the zero-range limit), and is consequently negligibly small[10]. The last 
term is the one we are seeking and was originally developed by Low[ll] in a limiting 
case for the deuteron, after the basic concept was sketched by Bohr [12]: 



AjE low = ™ 2 \uo)r J d , x J d3y w _ ^ Mf); B (f _ ^ x Jm + (g) 

A more convenient representation can be obtained by dividing both sides of this 
equation by the expression for the Fermi hyperfine energy given by Eqn. (1). Since 
the Wigner-Eckart theorem guarantees that Eqn. (8) must be proportional to a ■ S/S 
(which cancels in the ratio), we arrive at a simple expression for the leading-order 
nuclear contribution to the hfs, which is one of our primary results. 

AE^/E F = -2m e aS Low , (9) 
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where 

6how = ~i6/^ / d * x I d3y ^~^ m^)>((^ -y) x ( io ) 

and a nuclear expectation value is required of the z (or "3" ) component of this vector 
in the nuclear state with maximum azimuthal spin (i.e. S z = S). The intrinsic size 
of the nuclear corrections is given by (—2m e aR) = —38 ppm [i?/fm], where [i?/fm] 
is the value of the Low moment in Eqn. (10) in units of fm. 



3 Nuclear Matrix Elements 



In order to evaluate Eqn. (10) it is necessary to assume a form for the nuclear charge 
and current operators. We have agreed to ignore terms of relativistic order, and 
this eliminates all but the usual impulse-approximation (i.e., one-body) terms for 
the charge operator, which contains finite-size distributions for the protons and neu- 
trons. Although the latter contributions are small, they have never been included in 
previous calculations and we will gauge their importance by including them in our 
calculation. The nuclear current operator can be written in terms of the dominant 
spin-magnetization current, the convection current (motional current of charged par- 
ticles) and meson-exchange currents (MEC). Isoscalar MEC are of relativistic order [5] 
and we ignore them. Isovector MEC are larger, but don't contribute to the deuteron 
because it is an isoscalar nucleus. Isovector MEC will contribute to the trinucleon 
systems, where the effect of this current on the isovector magnetic moment is about 
the same size (15%) as our expansion parameter [13]. Because parts of these currents 
(in particular the isobar part) are difficult to treat quantitatively and because of their 
relative smallness, we will ignore MEC in calculating the Low moments in this initial 
effort to understand hfs in the trinucleon systems. 

Each of the charge and current operators that we use are therefore given by sums 
of one-body operators, and their resulting product in the Low-moment expression 
in Eqn. (10) can be written as a sum of one-body terms plus a sum of two-body 
terms. Using a transparent notation for this decomposition (5^ ow = 5^ in + 5^ w and 

^Low = ^i?ow + ^Lo^T) we fi 11 ^ that the one-body term is given for all nuclei by the 
spin-magnetization current in the form 

e - ms£ t (H*) < + m% t (^f) < ■ (») 

where 

(r>{5 = Jd 3 xJ d 3 y p p ch (x) p p M (y)\x - y\ = 1.086(12)fm (12a) 
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and 

(r)$ = Jd 3 xJ d 3 y p n ch {x) p n M {y)\x - y\ (126) 

determine the proton and neutron parts, respectively, of the one-body current. The 
quantities 7* and <7j are the (Pauli) isospin and spin operators for the ith nucleon, 
fP ch , p p M , and Pm are the proton charge and magnetic densities and the neutron mag- 
netic density (normalized to 1), while p™ h is the neutron charge density (normalized 
to 0). Note that the quantities (r)?^ and (r)" 2 " are the proton and neutron Zemach 
moments[8, 14], and we have listed in Eqn. (12a) the value of the proton Zemach mo- 
ment recently determined directly from the electron-scattering data for the proton[15] 
(the neutron has not yet been evaluated). In numerical work described below we 
will use simple forms for the neutron and proton form factors: a dipole form for 
the proton charge and magnetic form factors and the neutron magnetic form factor 
{Fu{q 2 ) = (1+^2/02)2 ) an d a modified Galster[16] form for the neutron charge form 

factor (Fq = nzf^mrp , where A = 0.0190 fm 2 ). To incorporate into our calculations 
the numerical value given by Eqn. (12a) we use f3 = 4.029 fm -1 , which reproduces 
that value for the dipole case. With this (3 the neutron Zemach moment has the value 
—0.042 fm, which is a very small correction to the proton. In Low's original work 
the nucleon Zemach moments were ignored (at that time no information existed that 
they were significant). 

The spin-isospin operators in Eqn. (11) are generators of SU(4) symmetry, and 
it is conventional [17] to decompose the wave functions of light nuclei with respect to 
that symmetry, which plays a significant role in understanding those systems. The 
dominant wave function component for the trinucleon systems (the S-state ~ 90%) 
is the product of a completely symmetric space wave function with a completely 
antisymmetric spin-isospin wave function. The next most important component is 
the D-state (~ 9%), and we will ignore the small remaining components for simplicity 
in the following discussion[10]. Treating only the proton term for the moment, we 

find the expectation value of J2f=i a 1 to ^ e ^z{l — for the deuteron, 

2S' Z (1 — \Pd) for the triton, and 2S Z ( — \Pd) for 3 He. The D-wave components tend 
to have the spin and orbital components anti-aligned, and this accounts for the sign 
of the Pd term. In 3 He the two protons tend to have their spins oppositely aligned, 
which accounts for the small 3 He Low moment. In the limit of exact SU(4) symmetry 
only the S-state contributes and the 3 He one-body term vanishes. 

There are three types of two-body Low moments: S-wave spin-magnetization 
terms, D-wave spin-magnetization terms, and (largely) D-wave convection current 
terms. For each of these types there are contributions from two protons, or from one 
proton and one neutron, or from two neutrons. We keep all terms but the convection 
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current contribution from two neutrons. The resulting nuclear operators are 

^Low = X! VjCij( x ij) ~ 7pij( X ij)(.3<7j • XijXij ~ &j) ) , (13a) 

Mat i-Li \ a / z 



where = x^- x (p^ — Pj), %ij = x% — Xj, x% is the coordinate of nucleon i, and p^ is 
the momentum of nucleon i. For simplicity we have not decomposed the radial (and 
isospin-dependent) functions C^, CV,-, and (7^- according to the types of nucleon that 
contribute. Explicit forms for these functions are given in Ref. [10]. In the limit of 
pointlike nucleons the radial part of each function becomes simply Xij. 



4 Results and Discussion 



The proton hfs has been discussed in detail recently[2, 15] and we have nothing 
additional to add. The recently evaluated proton Zemach moment was discussed in 
the text, and it leads to a —58.2(6) kHz contribution to the hydrogen hfs, which equals 
—41.0(5) ppm. When added to the usual QED and recoil corrections [1, 2, 15] there 
is a 3.2(5) ppm discrepancy with experiment, which can be attributed to hadronic 
polarization[18, 19] and (possibly) uncalculated recoil corrections. 

The deuterium, tritium, and 3 He + results were calculated using the (second- 
generation) AV18 potential [20], together with (for 3 H and 3 He) an additional TM' 
three-nucleon force[21] whose short-range cutoff parameter had been adjusted to pro- 
vide the correct binding energies. Individual one-body (labelled "nucleon" ) and two- 
body (labelled "Low") terms are tabulated together with their total in Table 2. The 
(approximate) SU(4) symmetry that dominates light nuclei [10, 17] provides an expla- 
nation for the relative sizes of the entries in this table, which we discuss next. Note 
that 3 He (which has proton number Z = 2) is uniformly enhanced by a factor of 
Z 3 = 8 contained in |0 n (O)| 2 in Eqn. (8). 

We expect (and verify) that interactions driven by the charge of the neutron will be 
significantly suppressed because the neutron is overall neutral. The neutron Zemach 
moment, for example, is ~ - 4% of that of the proton, and this greatly suppresses the 
neutron's contribution to the one-body term. The two protons in 3 He have their spins 
anti-aligned in the SU(4) limit, and this cancellation leads to a small net result for the 
one-body part. The protons in H and 3 H make comparable one-body contributions, 
since the proton in 3 H carries the entire spin in the SU(4) limit. 
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Table 2: Nuclear corrections to Is hyperfine structure in light hydrogenic atoms. 
"Nucleon" refers to the one-body part, and "Low" refers to the correlation or two- 
body term, and "Total" refers to their sum. All entries are given in kHz. 



H 


2 H 


3 H 


3 He+ 


Nucleon 


Nucleon Low Total 


Nucleon Low Total 


Nucleon Low Total 


-58.2(6) 


-41.1 87.3 46.2 


-50.6 -9.6 -60.1 


13.9 1428 1442 



The two-body terms that couple the neutron charge and the proton magnetic 
moment are suppressed for the reason discussed above. In addition the two-body 
convection current has no contribution from the dominant S-state and is therefore 
negligible for hydrogen isotopes, where one (charged) nucleon must be a neutron. 
Only for the two protons in 3 He is this interaction non-negligible, but still small. 
In 3 H the neutron spins are anti- aligned in the SU(4) limit, which suppresses the 
normally dominant proton charge - neutron magnetic moment contribution to about 
20% of the one-body part. In 3 He those terms add for the two protons, leading to 
an even larger result. Thus the qualitative features of the results in Table 2 can be 
understood in terms of (approximate) SU(4) symmetry. We quantify these qualitative 
observations in the following paragraph. 

The neutron's contribution to the nucleon one-body term is very small, except for 
3 He where it is ~ —40% of the very suppressed proton contribution. In deuterium the 
largest correction to the dominant proton charge - neutron magnetic moment scalar 
interaction (the first of the two terms in Eqn. (13a)) is the corresponding tensor 
interaction (the second term in Eqn. (13a)) and amounts to only 4%. In tritium the 
proton charge - neutron magnetic moment scalar interaction is suppressed by SU(4) 
symmetry, which enhances the relative contribution of the corresponding tensor term 
(~ 30%) and of the scalar neutron charge - proton magnetic moment term (~ 20%). 
The convection current contribution to both of these hydrogen isotopes is very small. 
The helium case is typified by many large contributions (in kHz) but the scalar proton 
charge - neutron magnetic moment term completely dominates. The scalar neutron 
charge - proton magnetic moment term is about 4% of the dominant interaction, the 
tensor terms are 1-2%, while the convection current of the two protons is about 5% 
of the dominant interaction, making it the largest of the corrections. 

Table 3 adds the results of Table 2 to the QED-only calculation for Is states, and 
expresses the differences with experiment as fractions of the Fermi energy. Results 
must be considered quite good, given the size of our hadronic expansion parameter. 
The deuterium case is particularly close to experiment, and this is likely due to the 
small binding energy, which tends to minimize relativistic corrections. The trinucleon 
cases range from very good in the 3 H case (~ 3% residue) to adequate in the 3 He case 
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(~ 20% residue). The large disparity in the two cases is undoubtedly due to missing 
MEC, particularly the isovector ones. Even this amount of missing strength is only 
slightly larger than our expansion parameter. 



Table 3: Difference between hyperfine experiments and hyperfine calculations for 
the Is state of light hydrogenic atoms, expressed as parts per million of the Fermi 
energy. The first line is the difference with respect to the QED calculations only, while 
the second line incorporates the hadronic corrections calculated above (the Zemach 
moment for hydrogen and nuclear corrections for the three nuclei). 



(£ms P - «)/M P pm) 



Theory 


H 2 H 3 H 3 He+ 


QED only 


-33 138 -38 -212 


QED + hadronic 


3.2(5) -3.1 1.2 -46 



Previous work on this topic is quite old[ll, 12, 22, 23, 24, 25], except for the 
deuterium [6] case. The older work relied on the Breit approximation for the electron 
physics, which is sufficient only for the leading-order corrections. It used an adiabatic 
treatment of the nuclear physics based on the Bohr picture of the nuclear hyperfine 
anomaly, which is far more complex than the treatment that we have presented. 
Uncalculated QED corrections and poorly known fundamental constants (such as a) 
led to estimates of nuclear effects that were many tens of ppm in error. Although 
the nuclear physics at that time was not adequate to perform more than qualitative 
treatments of the trinucleons, the SU(4) mechanism was known and this allowed 
a qualitative understanding. The only previous attempt to treat the three nuclei 
simultaneously was in Ref. [25]. They found nuclear corrections of about 200 ppm 
for deuterium, 20 ppm for 3 H, and —175 ppm for 3 He + . Except for the deuterium 
case (which involves significant cancellations) this has to regarded as quite successful, 
given the knowledge available at that time. 



5 Conclusions 



We have performed a calculation of the nuclear part of the hfs for 2 H, 3 H, and 3 He + , 
based on an expansion parameter adopted from x?T, a unified nuclear model, and 
modern second-generation nuclear forces. This is the first such calculation, and the 
results are quite good. Details of the results can be understood in terms of the 
approximate SU(4) symmetry that dominates the structure of light nuclei. 
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